Maximally Robust Unravelings of Quantum Master Equations 
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The stationary solution p of a quantum master equation can be represented as an ensemble of 
pure states in a continuous infinity of ways. An ensemble which is physically realizable through 
monitoring the system's environment we call an 'unraveling'. The survival probability S{t) of an 
unraveling is the average probability for each of its elements to be unchanged a time t after cessation 
of monitoring. The maximally robust unraveling is the one for which S{t) remains greater than the 
largest eigenvalue of p for the longest time. The optical parametric oscillator is a soluble example. 
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It is well-known that open quantum systems generally 
become entangled with their environment. Thus they 
do not remain pure, but become mixed. Nevertheless, 
pure states are used in the analysis of open quantum 
systems, in at least two distinct contexts. The first con- 
text is that of decoherence (see for example Refs. 
The pure states which are studied in this field are those 
which are relatively long-lived (or decohere slowly) un- 
der a particular irreversible evolution. In Rcf. m they 
are called stable states. The second context is that of 
continuous quantum measurement theory [^,^. Under 
some situations, the perfect monitoring of the system's 
environment collapses the system to a pure state which 
undergoes a stochastic evolution, called a quantum tra- 
jectory in Ref. Q. The ensemble of quantum trajectories 
reproduces the deterministic dissipative evolution, and so 
is called an unraveling of that evolution. 

In this letter we unite these two threads in the idea 
of a maximally robust unraveling. We consider systems 
obeying a quantum master equation, which is the most 
general generator of deterministic Markovian evolution 
and is used in many areas of physics. We further assume 
the system to have a unique stationary state matrix p^o ■ 
Monitoring the bath to which the system is coupled will 
give unravelings which are generated by a (nonlinear) 
stochastic Schrodinger equation (SSE). Different ways of 
monitoring the bath have different SSEs, but each SSE 
has a stationary ensemble of pure state solutions which 
represents Poo- Here we calculate the average survival 
probability S{t) of these physically realizable ensembles, 
rather than the survival probability of a single arbitrarily 
chosen pure state. As we explain below, the maximally 
robust unraveling is that for which S(t) remains higher 
than the largest eigenvalue of poo for the longest time. 
As an example we consider continuous Markovian unrav- 
elings of an optical parametric oscillator. 

Unraveling the master equation. — The most general 
form of the quantum master equation is j^J^ 



K 

p = -i[H, p] + V Ckpcl - \{clck,p} = Cp 



where {a, 5} = a6 + ba. The unique stationary state 
is defined by Cpoa = 0. This requires that C be time- 
independent, which implies that one can only remove a 
Hamiltonian term by moving into its interaction picture if 
the remainder of the superoperator C is time-independent 
in that interaction picture. 

It is now well-known that there are many (in fact con- 
tinuously many) different unravelings for a given master 
equation ||^ . If we represent the pure state of the system 
by the projector P, we can write any SSE unraveling the 
master equation (|^) as 



dP ^dt{C + U) P. 



(2) 



Here U (which depends on the operators Cfc) is a stochas- 
tic superoperator which is, in general, nonlinear in its 
operation on P. It is constrained by the following two 
equations which must hold for arbitrary projectors P 

[{C + U)P]P + dt[UP][UP] = (1- P)[{C+U)P] (3) 
E[UP] = 0. (4) 



The former property ensures that P 
jector and the latter that 

dE[P] = CE[P]dt, 



dP remains a pro- 



(5) 



where E denotes the ensemble average with respect to 
the stochasticity of U. This stochasticity is evidenced by 
the necessity of retaining the term dt[UP] [UP] in Eq. 

For simplicity we will call U an unraveling. Each un- 
raveling gives rise to an ensemble of pure states 



(6) 



where Pi are the possible pure states of the system at 
steady state, and Wi are their weights. For master equa- 
tions with a unique stationary state p^o, Wi is equal to 
the proportion of time the system spends in state Pi . The 
ensemble represents poo in that 



fc=i 



(1) 



(7) 
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as guaranteed by Eq. (||). Note that the projectors Pf^ 
wiU not in general be orthogonal. 

Robustness. — Let the system evolve under an unravel- 
ing U from an initial state at time — oo to the stationary 
ensemble at time 0. It will then be in the state P^^ with 
probability w'f . If we cease to monitor the system at 
this time then the state will no longer remain pure, but 
rather will relax toward poo under the evolution of Eq. |^). 
The robustness of a particular state Pf^ is measured by its 
survival probability SY{t). This is the probability that 
it would be found (by a hypothetical projective measure- 
ment) to still be in the state Pf at time t 

5f(t)=Tr[Ffcxp(/:i)^f]. (8) 

The average survival probability of the ensemble re- 
sulting from the unraveling lA is therefore 

S''{t)^Y.^Ys'Y^t). (9) 

In the limit t — > cxd the ensemble average survival prob- 
ability will tend towards the stationary value 6*^(00) = 
Tr[p^] which is independent of the unraveling U and is a 
measure of the mixedness of poo ■ An alternate measure 
of the mixednes of poo is its largest eigenvalue 

A= hm (Tr[p^])^/" = max{A : PooQa = AQa}. (10) 

n — *oo 

In the absence of any monitoring of the bath, the pro- 
jector Qa would be one's best guess for what pure state 
the system is in at steady state. The chance of this guess 
being correct is simply A and is obviously independent of 
time. Since S''^(0) = 1 > A > S'(oo), there will be some 
time t' at which the ensemble average survival proba- 
bility S^{t') becomes less than or equal to A. That is, 
given that we know by monitoring the environment the 
state Pf^ which is occupied at time then at time t' this 
state ceases (on average) to be any better than Qa as an 
estimate of which pure state is occupied by the system. 
In effect, the ensemble is obsolete at time t' , so it is 
natural to define the survival time of E^ as 

= miii{t : S'"(i) = A}. (11) 

Note that this expression depends on the entire superop- 
erator C, and so is sensitive to the interplay of Hamil- 
tonian and irreversible terms. This is as opposed to the 
reciprocal rate of decay of (t) at i = which depends 
only on the non-Hamiltonian terms in £ Q . 

The survival time quantifies the robustness of an 
unraveling U. Let the set of all unravelings be denoted 
J. Then the set of maximally robust unravelings Jm is 

jj^i = {ne J ■.t'^ >t'^ yu e J}. (12) 

Even if Jj\/ has many elements TZi, 7^2, • ■ these different 
unravelings may give the same ensemble E^ = E'^'^ — 
= . . .. There is clearly a strong reason to regard 



E^ as the most natural ensemble representation of the 
stationary solution of a given master equation. 

Continuous Markovian unravelings. — The formalism 
we have presented is well-defined but is difficult to apply 
because of the size of the set J of allowed unravelings. Al- 
though the stochasticity in the superoperators U can al- 
ways be written in terms of quantum jumps, these jumps 
range in size from being infinitesimal, to being so large 
that the system state after the jump is always orthogo- 
nal to that before the jump Also the unravelling need 
not be Markovian, even though the master equation is [Q. 
For this reason it is useful to consider a smaller (but still 
continuously infinite) set J' containing only continuous 
Markovian unravelings. A continuous (but not differen- 
tiable) time evolution arises from infinitely small (and 
infinitely frequent) jumps. In this case the probability 
distribution for the pure states obeying the SSE satis- 
fies a Fokker-Planck equation These unravelings are 
thus the natural ones to consider for quantum systems 
expected to show quasi-classical behaviour ||] . 

The elements W of J' can be written as 

K 

U'dt = J2'H[dWk{t)ck]. (13) 
fc=i 

Here Ti.[r] is a nonlinear superoperator defined (for an 
arbitrary operator r) by 

'H[r]p = rp + pr^ — Tr[rp + pr^]p (14) 

and the dWk (t) are the infinitesimal increments of a com- 
plex multi-dimensional Wiener process [|j satisfying 

dWj it)dW^ (t) = dtSjk , dWj it)dWk {t) = dtujk , (15) 

where the Ujk are arbitrary complex numbers obeying 
\ujk\ < 1 and Ujk = Ukj. 

An example. — As an application for our concept of 
robust unravelings we consider the master equation 

Cp - -(x/4)[a' - (at)2,p] + apa^ - \{a^a,p}. (16) 

Here a is the annihilation operator for a harmonic os- 
cillator, satisfying [a, a^] = 1. For |x| < 1 this mas- 
ter equation is a good model for a degenerate optical 
parametric oscillator below threshold [Q. Time is be- 
ing measured in units of the cavity lifetime. The free 
Hamiltonian Hq = Loa'^a is omitted because we are work- 
ing in its interaction picture. This is necessary for the 
Hamiltonian term in Eq. (|l6|), which arises from driving 
a x^^^ crystal with a classical field at frequency 2ijj^ to 
be time-independent. The damping term remains time- 
independent in the interaction picture. 

For Ixl < 1 the equation Cpoo = has a solution, and 
it is a Gaussian state. That is to say, the symmetrized 
moments for the quadrature operators x = a -\- a'^ and 
y = —ia + ia^ are those of a bivariate Gaussian distri- 
bution (which is in fact the Wigner function for the 
state). Such states are characterized by a mean vector 
and a covariance matrix which we denote 
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M = 



{x^) \{xy + yx) \ _ 



^{xy + yx) 



7 f3 
(i a 



(17) 
(18) 



For all Gaussian quantum states the second-order mo- 
ments obey the following inequality 



(Tr[p2]) 2 = dct(M) = a7 - > 1, 



(19) 



where the equality holds only for pure states. For we 
find Moo = cliag[(l + (1 - x)'^] and Moo = 0. The 

maximum eigenvalue for this state is 



A = 2 



l + (l-x')-l/2 



(20) 



The continuous Markovian SSEs which unravel the 
master equation ( [T^ ) are of the form 



dP = CPdt + rL[dW{t)a]P, 



(21) 



where dW*{t)dW{t) ~ dt. Because there is only one ir- 
reversible term in the master equation, these unravelings 
are parameterized by the single complex number u satis- 
fying |u| < 1 and defined by 



dW{t)^ = udt ^ {r + ih)dt. 



(22) 



The different unravelings (hereafter denoted by u) could 
be physically realized by dividing the cavity output into 
two beams so that two different quadratures may be mea- 
sured simultaneously by homodyne detection. 



To find the ensemble E" we need to solve the SSE (|2lD . 
This can be done by assuming a Gaussian (as defined 
above) initial state. Then, following the method of Ref. 
pO| , we can show that this state will remain Gaussian, 
with its moments obeying the differential equations 

dx = -(1 + x){x/2)dt + Re {dW{t){-f - 1 ~ ioi)\ , (23) 

dy = -(1 - x)(y/2)dt -I- Re {dW(t){li - ia + i)} , (24) 
7 = -(1 + X)7 + 1 + (1/2) {-(1 + r)(7 - 1)' 

- (1 - r)/32 + 2/i(7 - l)/3} , (25) 

d = -(1 - X)« + 1 + (1/2) {-(1 - r){a - If 

- {l+r)f3^ + 2h{a~l)/3} , (26) 
$ = -p+ (1/2) {-(1 + r)(7 - 1)13 - (1 - r)ia - l)/3 

+ h[p' + - l)ia - 1)]} . (27) 

Here the terms arising from the unraveling U = 
'H[dW{t)a] have been enclosed in curly brackets. 

What is remarkable about these equations is that the 
second-order moments evolve deterministically. It can be 



verified numerically that they have a stationary solution 
Mq which satisfies the equality in Eq. (|l^), as expected. 
For a given unraveling u ^ r + ih, the system will evolve 
from an arbitrary state at time — oo into a unique en- 
semble at time consisting of Gaussian pure states PS, 
all with the same second-order moments Mg but distin- 
guished by their first-order moments fl. Thevector jl thus 
acts as the index i for the ensemble as in Eq. (^). Fur- 
thermore, because Mq is constant, Eqs. ( [2^ ) and ( p4[ ) for 
the first-order moments are those of a two-dimensional 
Ornstein-Uhlenbeck equation The stationary prob- 
ability distribution wp for fl is therefore Gaussian. By 
virtue of Eq. (^ the second-order moments of this distri- 
bution plus those of the pure states Pp must equal those 
of Poo- That is, is a bivariate Gaussian in fl with 
mean and covariance matrix Moo — Mq . 

Due to the nonlinearity of the equations of motion for 
M", there is no analytical expression for the stationary 
solutions Mq in terms of u. However the inverse prob- 
lem is easy to solve because the three moment equations 
(p5|)-(p7[) with the left-hand sides equal to zero are linear 
in r and h. That is, given the moments Mq satisfying 
det(MQ) = 1, these three equations are consistent and 
yield values of r and h in terms of 7, a, Of course 
only those solutions satisfying + < 1 correspond to 
physically realizable unravelings. This easily yields the 
equation for the set of allowed moments Mq . This region 
is plotted in Fig. 1 for x = 0.9. Only two axes (namely 
/3o,7o) are needed because the equation det(MQ) = 1 
implies ajf = mf + 1]/^^. 

This plot of the physically realizable ensembles con- 
tains information about the system not revealed by pre- 
vious work. In particular the set of these ensembles is 
smaller than the set of all ensembles of Gaussian pure 
states with fixed second-order moments but varying first- 
order moments which reproduce the stationary state poo • 
The latter set of second-order moments M is defined by 
the two equations det(M) = 1 and det(Moo - M) > 0. 
This region is also plotted in Fig. 1. 

To find the survival probability we need the time- 
evolved state e^^Pp, which is still Gaussian. From 
Eqs. (p3|)-(p7|) with the terms in curly brackets removed 
its moments are found to be 



Pt^Vtfi, MJ" - VtMJJVf + M00- V^MooVf, (28) 

where Vj — diag(w-|-, f_), where v± — exp ( — (l±x)t/2). 
Using the above results and the properties of Wigner 
functions M we find 



d'flw-pTT[Pi:exp{£t)Pi:] = {det(i[(l-V,)(Moo-M«)(l-Vn + Mr + Mff])} 



^1/2 



(29) 



(l + x)7o" 



1-X 



X 



(00 



i«\2 



(1 + xH 



-1/2 



(30) 
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It is easy to prove that for < x < 1 this expression 
is maximized at all values of < f < oo by choosing 
7o = (1 + x)~^ (its maximum allowed value) and Pq = 0. 
From Eqs. (|2^)~(^) this is achieved only for the unrav- 
eling u — —1. For this maximally robust unraveling the 
survival probability simplifies to 



^2g-(l-x)t/2 ■ 



(31) 



From Eqs. (|ll|) and ( |20| ) the maximum ensemble average 
survival time is thus 
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4x' ( 2 + X' 



(32) 



This is plotted as a function of x in Fig- 2, along with 
other parameters of interest. Note that the survival time 
grows with x from = 2 In 2 at x = and diverges 
as (1 — x)~^ as X ^ 1- This is characteristic of criti- 
cal slowing down. The increase in with x does not 
mean that the ensemble becomes more robust as x 
increases; on the contrary, for fixed t, S^{t) decreases 
monotonically with x- For large x the survival probabil- 
ity S'^{t) decays more quickly, but remains greater than 
the maximum eigenvalue A (which decreases rapidly for 
large x) for a longer time. 

In Conclusion, maximally robust unravelings are a 
quantitative tool for investigating open quantum sys- 
tems, and give insights quite distinct from those offered 
by traditional techniques. Our work differs from previ- 
ous work in that we have considered the average survival 
probability S(t) for an ensemble which is physically real- 
izable from monitoring the environment, rather than the 
survival probability (or change in entropy) of an arbitrar- 
ily chosen single pure state. Also our definition of the 
survival time r as the time at which the ensemble result- 
ing from monitoring the environment becomes obsolete is 
fundamentally motivated, unlike other definition such as 
= As{t)\t=o or S{t) = 1/e. We have shown that 
for a particular class of unravelings, those which are con- 
tinuous and Markovian, an analytical solution is possible 
for a simple system. Our technique is easily generalizable 
to any system which can be linearized, although a partly 
numerical solution may be necessary. 



FIG. 1. Plot of the allowed values of the second-order mo- 
ments 7 = (a;^) and f3 = {xy + yx) /2 held by all members 
of an ensemble of Gaussian pure states which represent the 
stationary state poo- The solid line is the boundary of the al- 
lowed region for ensembles which are physically realizable 
through Eq. (pl|), while the dashed line is for unconstrained 
ensembles. The star indicates the position of the most robust 
ensemble E^. The threshold parameter x equals 0.9. 
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FIG. 2. Plot of (a) the average survival time (in units of 
the decay time) of the maximally robust ensemble E^; (b) the 
y- variance Ooo for the stationary state poo; (c) the j/- variance 
Of J' for the members of E^; (d) the largest eigenvalue A of 
Poo; and (e) 5(00) — Tr[p^] versus threshold parameter x- 
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